clc;
syms n k positive;
syms p real;
assume(n,'integer');
disp("----------------------下面是判定方法（3）--------------------------------------------")
a=(symprod(2*k-1,k,1,n)/symprod(2*k,k,1,n))^p*((x-1)/2)^n;
F=simplify(subs(a,n,n+1)/a),
L=simplify(limit(F,n,inf))%书本方法（3）判定
disp("----------------------下面是判定方法（4），且x=2-------------------------------------")
b=(symprod(2*k-1,k,1,n)/symprod(2*k,k,1,n))^p;
L=limit(n*(b/subs(b,n,n+1)-1),n,inf)